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Metal-insulator transition in a one-dimensional magnetic-non-magnetic superlattice 
geometry: Existence of multiple mobility edges 
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Based on Green's function formalism, the existence of multiple mobility edges in a one-dimensional 
magnetic-non-magnetic superlattice geometry in presence of external electric field is predicted, and, 
it leads to the possibility of getting a metal-insulator transition at multiple values of Fermi energy. 
The role of electric field on electron localization is discussed for different arrangements of magnetic 
and non-magnetic atomic sites in the chain. We also analyze that the model quantum system can 
be used as a perfect spin filter for a wide range of energy. 
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I. INTRODUCTION 

Quantum transport in low-dimensional systems has 
been a topic of interest within the past few decades due to 
its potential applicability in the field of nanoscience and 
nanotechnology. Exploitation of the spin degree of free- 
dom adds a possibility of integrating memory and logic 
into a single device, leading to remarkable development 
in the fields on magnetic data storage application, de- 
vice processing technique, quantum computation^, etc. 
Naturally a lot of attention has been paid to study spin 
transport in low-dimensional systems both from experi- 
mental^r^ and theoretical^"— points of view. 

The understanding of electronic localization in low- 
dimensional model quantum systems is always an inter- 
esting issue. Whereas, it is a well established fact that 
in an infinite one-dimensional (ID) system with random 
site potentials all energy eigenstates are exponentially lo- 
calized irrespective of the strength of randomness due to 
Anderson localizationi^, there exists another kind of lo- 
calization, known as Wannier-Stark localization, which 
results from a static bias applied to a regular ID lat- 
tice, even in absence of any disorder—. Till date a large 
number of works have been done to explore the under- 
standing of Anderson localization and scaling hypoth- 
esis in one- and two-dimensional systemsi^. Similarly, 
Wannier-Stark localization has also drawn the attention 
of many theoristsi^Ti^^ as well as experimentalists^i. For 
both these two cases, viz, infinite ID materials with ran- 
dom site energies and ID systems subjected to an exter- 
nal electric field, one never encounters any mobility edge 
i.e., energy eigenvalues separating localized states from 
the extended ones, since all the eigenstates are localized. 
But there exist some special types of ID materials, like 
quasi-periodic Aubry- Andre model and correlated disor- 
dered systems where mobility edge phenomenon at some 
particular energies is obtained^^"— . Although the studies 
involving mobility edge phenomenon in low-dimensional 
systems have already generated a wealth of literature^"— 
there is still need to look deeper into the problem to ad- 



dress several interesting issues those have not yet been 
explored. For example, whether the mobility edges can 
be observed in some other simple ID materials or the 
number of mobility edges separating the extended and 
localized regions in the full energy band of an ID mate- 
rial can be regulated, are still to be investigated. 

To address these issues in the present article we in- 
vestigate two-terminal spin dependent transport in a 
ID mesoscopic chain composed of magnetic and non- 
magnetic atomic sites in presence of external electric 
field. To the best of our knowledge, no rigorous effort 
has been made so far to explore the effect of an ex- 
ternal electric field on electron transport in such a fD 
magnctic-non-magnctic superlattice geometry. Here we 
show that, depending on the unit cell configuration, a 
ID superlattice structure subjected to an external elec- 
tric field exhibits multiple mobility edges at different val- 
ues of the carrier energy. We use a simple tight-binding 
(TB) framework to illustrate the model quantum sys- 
tem and numerically evaluate two-terminal spin depen- 
dent transmission probabilities through the superlattice 
geometry based on the Green's function formalism. From 
our exact numerical analysis we establish that a sharp 
crossover from a completely opaque to a fully or partly 
transmitting zone takes place which leads to a possibility 
of tuning the electron transport by gating the transmis- 
sion zone. In addition to this behavior we also show that 
the magnetic-non-magnetic superlattice structure can be 
used as a pure spin filter for a wide range of energy. These 
phenomena enhance the prospect of such simple superlat- 
tice structures as switching devices at multiple energies 
as well as spin filter devices, the design of which has sig- 
nificant impact in the present age of nanotechnology. 

With an introduction in Section I, we organize the pa- 
per as follows. In Section II, first we present the model, 
then describe the theoretical formulation which include 
the Hamiltonian and the formulation for transmission 
probabilities through the model quantum system. The 
numerical results are illustrated in Section III and finally, 
in Section IV, we draw our conclusions. 



II. THEORETICAL FRAMEWORK 

Let us start with Fig. [1] where a ID mesoscopic chain 
composed of magnetic and non-magnetic atomic sites 
is attached to two semi-infinite ID non-magnetic elec- 
trodes, namely, source and drain. The chain consists 
of p {p being an integer) number of unit cells in which 
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FIG. 1: (Color online). A ID mesoscopic chain composed 
of magnetic (filled magenta circle) and non-magnetic (filled 
green circle) atomic sites is attached to two semi-infinite ID 
non-magnetic metallic electrodes, namely, source and drain. 

each unit cell contains n and m numbers of magnetic 
and non-magnetic atoms, respectively. Both the chain 
and side-attached electrodes are described by simple TB 
framework within nearest-neighbor hopping approxima- 
tion. 

The Hamiltonian for the entire system can be written 
as a sum of three terms as. 



H = H, + Hi+ Ht, 



(1) 



The first term represents the Hamiltonian for the chain 
and it reads 
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Here, e,; refers to the on-site energy of an electron 
at the site i with spin a (t,i), t is the nearest- neighbor 
hopping strength, Cj^-t {cia) is the creation (annihilation) 
operator of an electron at the ith site with spin a and hi 
is the strength of local magnetic moment where hi = 
for non-magnetic sites. The term hi. a corresponds 
to the interaction of the spin of the injected electron 
with the local magnetic moment placed at the site i. 
The direction of magnetization in each magnetic site is 
chosen to be arbitrary and specified by angles 0i and 
(j)i in spherical polar co-ordinate system for the ith 
atomic site. Here, 9i represents the angle between the 
direction of magnetization and the chosen Z axis, and 4>i 
represents the azimuthal angle made by the projection 
of the local moment on X-Y plane with the X axis. 
In presence of bias voltage V between the source and 



drain an electric field is developed, and therefore, the 
site energies of the chain becomes voltage dependent. 
Mathematically we can express it as e^ = e^* + ^i{V), 
where e^ is the voltage independent term. The voltage 
dependence oi ei{V) reflects the bare electric field in the 
bias junction as well as screening due to longer range 
electron-electron interaction. In the absence of such 
screening the electric field varies uniformly along the 
chain and it reads ei{V) = V/2 — iV/{N + 1), where 
N corresponds to the total number of atomic sites in 
the chain. In our present work, we consider both the 
linear and screened electric field profiles. As illustrative 




FIG. 2: (Color online). Voltage dependent site energies in 
a ID chain considering 100 atomic sites for three diflterent 
electrostatic potential profiles when the bias voltage V is set 
equal to 1. 

example, in Fig. [5] we show the variation of voltage 
dependent site energies for three different electrostatic 
potential profiles for a chain considering 100 atomic sites 
and describe the nature of electronic localization for 
these profiles in the forthcoming section. 

The second and third terms of Eq. [1] describe the TB 
Hamiltonians for the ID semi-infinite non-magnetic elec- 
trodes and the chain-to-electrode coupling. These Hamil- 
tonians are written as follows. 



H, 



and. 
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^ cleiCn + ^ [cltiCn+1 + h.c] 



(3) 



Htun — Htun,S + Htun.D 

= Ts[c\cQ^h.c.]+Td[c^j^CN+i + h.c]. (4) 

The summation over S and D in Eq. [3] implies the incor- 
poration of both the two electrodes, viz, source and drain. 
£; and ti stand for the site energy and nearest-neighbor 
coupling, respectively. The electrodes are directly cou- 
pled to the chain through the lattice sites 1 and TV, and 
the coupling strengths between these electrodes with the 
chain are described by t^ and r^, respectively. 

To obtain spin resolved transmission probabilities of an 
electron through the source-chain-drain bridge system, 
we use Green's function formalism. The single particle 
Green's function operator representing the entire system 
for an electron with energy E is defined as. 



G= {E-H + ir,y^ 



(5) 



where, ry — >■ 0+. 

Following the matrix form of H and G the problem 
of finding G in the full Hilbert space H can be mapped 
exactly to a Green's fmiction G^^^ corresponding to an 
effective Hamiltonian in the reduced Hilbert space of the 
chain itself and we have, 



■^eff 
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where, 



flipping can be explained from the following arguments. 
The operators (7+ (= a^+icy) and a- (= a^ —io'y) asso- 
ciated with the term hi.a in the TB Hamiltonian Eq. [5] 
are responsible for the spin flipping, where ff being the 
Pauli spin vector with components a^ , Ty and cr^ for the 
injecting electron. In our present model since we consider 
that all the magnetic moments are aligned along +Z di- 
rection, the term hi. a (= hi^ax + hiyGy -\- hiz<Jz) gets 
the form hizGz, and accordingly, the Hamiltonian docs 
not contain Ux and Cy and so cr+ and (t_ do not appear. 
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These Sg and So are the self-energies introduced to in- 
corporate the effect of coupling of the chain to the source 
and drain, respectively. Using Dyson equation the ana- 
lytic form of the self energies can be evaluated as follows. 
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where, Ci = iE-ei)/2 



i^ 



{E - eiY/A. 



Following Fisher-Lee relation, the transmission proba- 
bility of an electron from the source to drain is given by 
the expression. 



T,,,=Tv[T''sg'Tig'']. 



(9) 



where, r^,^s 's are the coupling matrices representing the 
coupling between the chain and the electrodes and they 
are defined as. 
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Here. 



S^ and S^' are the retarded and advanced self- 



energies associated with the fc-th (fc = 5, D) electrode, 
respectively. 

Finally, we determine the average density of states 
(ADOS), p{E), from the following relation. 



P{E) 



1 



Im[Tr[a] 



(11) 



In what follows we limit ourselves to absolute zero tem- 
perature and use the units where c — e^h~l. For the 
numerical calculations we set t = 1, e^ = OVi, /i; = 1 
for the magnetic sites, 9i = (j)i = 0, e; = 0, i; = 1 and 
Ts = Td =0.8. The energy scale is measured in unit of t. 

III. NUMERICAL RESULTS AND DISCUSSION 

Throughout our numerical calculations we assume that 
the magnetic moments are aligned along -\-Z direction 
{6i = (pi =0), which yields vanishing spin flip transmis- 
sion probability, viz, T^^_ = T^^ = 0, across the bridge sys- 
tem. The net transmission probability is therefore a sum 
T{E) = T-^-^{E)+Tii{E), and the origin of this zero spin 
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FIG. 3: (Color online). Transmission probability T and 
ADOS as a function of energy for a ID magnetic-non- 
magnetic superlattice geometry considering a linear bias drop 
along the chain, as shown by the pink curve in Fig. (2] where 
(a)-(c) correspond to the results for three different values of 
bias voltage V. 

which leads to the vanishing spin flip transmission proba- 
bility across the ID chain. Below, we address the central 
results of our study i.e, the possibility of getting multiple 
mobility edges in ID magnctic-non-magnetic superlattice 
geometries and how such a simple model quantum sys- 
tem can be used as a perfect spin fllter for a wide range 
of energy. 

In Fig. [3] we show the variation of total transmission 
probability T along with the average density of states 
for a ID magnetic-non-magnetic superlattice geometry 
considering a linear bias drop. Here we consider a 400- 
site chain in which each unit cell contains one magnetic 



and four non-magnetic sites and the results are shown for 
three different bias vohages. For the particular case when 
the chain is free from external electric field i.e., V = 
electronic conduction through the bridge takes place for 
the entire energy band as shown in Fig. ^a) which pre- 
dicts that all the energy eigenstates are extended in na- 
ture. The situation becomes really very interesting when 
the superlattice geometry is subjected to an external elec- 
tric field. It is illustrated in Figs. [S^b) andlSjc). From 
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FIG. 4: (Color online). Transmission probability T and 
ADOS as a function of energy for a ID magnetic-non- 
magnetic superlattice geometry when the electrostatic poten- 
tial profile varies following the green curve shown in Fig. (2] 
where (a)-(c) represent the identical meaning as in Fig. [31 

these spectra we notice that there are some energy re- 
gions for which the transmission probability completely 
drops to zero which reveals that the eigenstates associ- 
ated with these energies are localized, and they are sepa- 
rated from the extended energy eigenstates. Thus, sharp 
mobility edges are obtained in the spectrum, and, the 
total number of such mobility edges separating the ex- 
tended and localized regions in a superlattice geometry 
in presence electric field strongly depends on the unit 
cell configuration and it can be regulated by adjusting 
the number of magnetic and non-magnetic sites. This 
phenomenon describes the existence of multiple mobility 
edges in a superlattice geometry under finite bias con- 
dition. Now if the Fermi energy is fixed at a suitable 



energy zone where T drops to zero an insulating phase 
will appear, while for the other case, where T is finite, a 
metallic phase is observed and it leads to the possibility 
of controlling the electronic transmission by gating the 
transmission zone. The width of the localized regions 
between the band of extended regions increases with the 
strength of the electric field as clearly shown by com- 
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FIG. 5: (Color online). Transmission probability T and 
ADOS as a function of energy for a ID magnetic-non- 
magnetic superlattice geometry when the electrostatic poten- 
tial profile varies following the blue curve shown in Fig. [51 
where (a)-(c) represent the identical meaning as in Fig. [S] 

paring the spectra given in Figs. [3l[b) andlUJe), and, for 
strong enough field strength almost all energy eigenstates 
are localized. In that particular limit metal-to-insulator 
transition will no longer be observed. 

The above results are analyzed for a particular (linear) 
variation of electric field along the chain. To explore the 
sensitivity of getting metal-to-insulator transition on the 
distribution of electric field, in Figs. S] and [5] we present 
the results for two different screened electric field pro- 
files taking the identical chain length. From the spectra 
we clearly observe that the width of the localized region 
gradually disappears with the flatness of the electric field 
profile in the interior of the bridge system. If the poten- 
tial drop takes place only at the chain-to-electrode inter- 
faces, i.e., when the potential profile becomes almost fiat 
along the chain the width of the localized region almost 



vanishes and the metal-to-insulator transition is not ob- 
served, as is the case for the zero bias limit. 

Finally, we illustrate how such a simple magnetic-non- 
magnetic superlattice geometry can be utilized as a per- 
fect spin filter for a wide range of energy in absence of any 
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FIG. 6; (Color online). T^t, T^ and ADOS as a function of 
energy for a ID magnetic-non-magnetic superlattice geometry 
in absence of external electric field. 

external electric field. As illustrative example, in Fig. [S] 
we present the transmission probabilities for up and down 
spin electrons together with the average density of states 
as a function of energy for a ID magnetic-non-magnetic 
superlattice geometry. From the spectra we observe that 
the up and down spin electrons follow two different chan- 
nels while traversing through the superlattice geometry, 
since the spin flipping is completely blocked for this con- 
figuration. This splitting of up and down spin conduc- 
tion channels is responsible for spin filtering action and 
the total number of these channels strongly depends on 
the unit cell configuration. From Figs. [6l^a) and (b) we 
clearly see that for a wide range of energy for which 



the transmission probability of up spin electrons drops 
to zero value, shows non-zero transmission probability of 
down spin electrons. Therefore, setting the Fermi energy 
to a suitable energy region we can control the transmis- 
sion characteristics of up and down spin electrons, and, a 
spin selective transmission is thus obtained through the 
bridge system. Before we end, we would like to point 
out that since the overlap between the up and down spin 
conduction channels depends on the magnitudes of the 
local magnetic moments, we can regulate the spin degree 
of polarization (DOP) simply by tuning the strength of 
these magnetic moments and for a wide range of energies 
it (DOP) almost reaches to 100%. Thus, our proposed 
magnetic-non-magnetic superlattice geometry is a very 
good example for designing a spin filter. 



A. Conclusion 

To conclude, in the present work we investigate in de- 
tail the spin dependent transport under finite bias con- 
dition through a ID magnctic-non-magnetic superlattice 
geometry using Green's function formalism. Wc use a 
simple TB framework to describe the model quantum 
system where all the calculations arc done numerically. 
From our exact numerical analysis we predict that in 
such a simple ID magnetic-non-magnetic superlattice ge- 
ometry multiple mobility edges separating the localized 
and extended regions are obtained in presence of exter- 
nal electric field and the total number of mobility edges 
in the full energy spectrum can be controlled by arrang- 
ing the unit cell configuration. This phenomenon reveals 
that the superlattice geometry can be used as a switching 
device for multiple values of Fermi energy. The sensitiv- 
ity of metal-to-insulator transition and vice versa on the 
electrostatic potential profile is thoroughly discussed. Fi- 
nally, we analyze how such a superlattice geometry can 
be utilized in designing a tailor made spin filter device 
for wide range of energies. Setting the Fermi energy at 
a suitable energy zone, a spin selective transmission is 
obtained through the bridge system. All these predicted 
results may be utilized in fabricating spin based nano 
electronic devices. 

The results presented in this communication are 
worked out for absolute zero temperature. However, 
they should remain valid even in a certain range of finite 
temperatures (~ 300 K). This is because the broaden- 
ing of the energy levels of the chain due to the chain-to- 
electrode coupling is, in general, much larger than that 
of the thermal broadening^. 
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